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Abstract. From the inception of radio astronomy, the study of the Interstellar Medium 
has been both aided and frustrated by one fact: we are right within it! Our embedded 
perspective is favourable to observing kinematic and structural details that are likely to 
remain inaccessible in external galaxies for some time, but not so to gaining a birds-eye 
view of our Galaxy's structure and motions on the largest scales. The Canadian Galac- 
tic Plane Survey (CGPS) is the original survey with the ability to image both details 
and the "big picture". We will broadly review what is known of the Milky Way, and 
focus on large-scale ISM structure and dynamics that the CGPS and immediate family 
of surveys depicts particularly well (e.g. spiral structure, the rotation curve, density 
waves, rolling motions, the warp & scalloping). We also highlight areas where puzzles 
still exist (e.g. outer spiral structure, the question of shocks and rolling motions in the 

Milky Way), and offer some new insights (e.g. multiple shocks in the 2 nc ' quadrant; a 
radially varying spiral pattern speed in the disc) that demonstrate what is possible with 
current and future high-resolution 21 cm surveys. 



1. Introduction and a Brief History 

Although the meeting that generated these proce edings coincides w ith the end of ob- 
serving for the Canadian Galactic Plane Survey dTaylor et al.ll2003h and celebrates the 
completion of radio astronomy's original "small-structure: big-picture" survey, it is 
more properly viewed as a rest-stop along the journey that all ISM astronomers share 
in. Discoveries made at the survey's end will undoubtedly be as profound as the ones 
we have made along the way and will continue to make. The journey for those who 
have been involved with the CGPS has also been marked with images of an extraordi- 
nary beauty only a well-sampled u, v plane could provide. This review will take stock 
of discoveries we have accumulated up to this point and synthesize a contemporary 
picture of the Milky Way's large-scale structure and dynamics. It follows a talk on the 
same, online at www.novell.brandonu.ca/~fostert/foster.ppt. As the author is a progeny 
of the CGPS, the review is naturally focussed on aspects of the Milky Way (hereafter 
MW) that are depicted well in it and its siblings VLA-GPS and Southern-GPS, which 
came together in 2001 to form the International Galactic Plane Survey, or IGPS. 

One theme that the reader may expect to encounter in pursuing or reviewing Galac- 
tic ISM structure and dynamics studies is an often frustrating interdependency between 
them. The cause is in the brightness temperature unit TB(v)dv (K), which descends 
from both line-of-sight (hereafter LOS) density variations «///(r) and inflections in the 
velocity fi eld dv/dr. Th is "degeneracy" frustrated the 1970's pioneers of the field; for 
tonjlT 



example. iBurtonl ( 197 ll) attempted to say something about the distribution of Galactic 
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hydrogen by first assuming its motion (rotation only), whereas Rohlfs' 1 19741 investiga- 
tion of the velocity field could only proceed by assuming some density structure for the 
disk. Along with the confusion from the overwhelming amount of emission, the nearly 
indistinguishable observational signatures of different density wave theories, and the 
lack of computing power at the time, this led to early efforts to model spiral structure 
with H i being largely abandoned by the late 1970's. The demise was foreshadowed by 
Frank Kerr's comments on Galactic H i spir al structure at the 1970 "Spiral Workshop" 
(an account of which is by ISimonsonlll970T) . where he introduced 21 cm emission as 
"the way to study the whole Galaxy", onl y to concede late r that "there is so much 21 cm 
emission that we're confused by it" (see iLockmanl 120021) . As you will see, the picture 
remains elusive, and to this day there are doubters that H i is at all a reliable tracer of 
spiral structure. 



2. Fundamental Constants 

The fundamental equation of Galactic structure relates the observed Doppler velocity 
along the line-of-sight vios to the Galactic orbital motion Q of an ISM element: 

vlos = Ro (R) - Q ] sin£ cosZ? (1) 

where the fundamental structural constant Ro is the Solar Galactocentric radius and the 
fundamental dynamical constant Qo = vq/Ro is the angular velocity of the Sun (with 
respect to the LSR). Mapping out 21 cm Galactic structure and dynamics requires each 
of Q and R to be measured independently for the element. Since vios is measured, Q 
is an observable quantity independent of R, provided both Qo and Ro are known. Only 
if perturbations in R and z are negligible is vlos then a direct observation of the true 
rotational frequency Q of the element with respect to the Sun and Eqn. [T]is independent 
of R. Unfortunately, these ideal conditions do not often happen in the plane. 

The angular velocity in the Solar neighborhood Qo is most often measured by 
observing kinematics of stellar tracers such as OB stars, Cepheids, open star clus- 
ters, or old red giants. All methods either rely on observables like the Oort con- 
stants, (Qo - A - B) or require the distance to each tracer (and hence rely on the 
constant Rq). It appears that differences in measured values for both constants are 
subject to the observer's bias more so than to random differences from the techniques 
used to measure them. Hence neither Qo or Rq have converged towards single val- 
ues with time. A gradual increase in both is shown in Figure \T\ which plots 32 pub- 
lished measurements of Ro and 23 for Qo- The data (Table 1, available online at 
www.novell.brandonu.ca/~fostert/tablel.pdf) are grouped together into one-year bins 
by calculating the uncertainty-weighted mean for each year (labelled in figure). The 
uncertainty for each bin is the mean published uncertainty that year divided by the 
square root of the number of measurements that year (in brackets on figure). The trend 
in Qo (Pearson coefficient r=0.87) is significant and likely attributable in part to the 
"bandwagon" effect where new values will tend to agree with previously published 
ones regardless of the method used to obtain them. A correlation of Ro with publi- 
cation year is also observe d (r-0.59; 95% level). This suggests that the same effect 
in Ro (previously noted by Reid1 ll993h has carried forward beyond 1993, although all 
data bins in Fig. Q] are within error of one another. An error- weighted mean for each 
constant Ro =8.0±0.4 kpc and Qo —28.7+ 1. 1 km s~ l kpc" 1 is calculated, though it is 
likely that neither sample is complete. For this review we adopt the modern constants 
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of ReidetaL | (|2009b) based on trigonometric parallax observations of star formation 
regions: Rq =8.4±0.6 kpc and Oo =30.3±0.9 km s" 1 kpc' 1 (or a Solar linear velocity 
of v -25410ns- 1 ). 




Figure 1. Published Galactic constants Rq («; left axis) and Q.q (□; right) over time. 



3. Galactic Dynamics 



3.1. Milky Way Circular Rotation 

The largest scale dynamics in the Milky Way are those that reflect the response of 
the matter to the overall gravitational potential Vo (R) created by the mass distribution 
(stars plus ISM). For circular rotation the potential is related to observables Q. and R 
by dVo/dR - Q. 2 R. Circular orbits in galaxies are found with fl oc R~ l or a. constant 
circular speed, and never at the extremes CI oc R~ 3 / 2 (Keplerian) and £1 - const (solid- 
body rotation). The actual shape of MW rotation near the Sun can be assessed with the 
Oort constants, or equivalently the frequency of radial oscillations about the circular 
orbit's guiding centre kq = ^-4B(A - B). The ratio of epicyclic-to-orbital frequency, 
k/Q is V2 for flat circular speed cur ves. This is in fact w hat is observed near the Sun: 



for example using open star clusters iLepine et al.1 (12008) find kq =43±5 km s kpc 



Then 43/30.3 ~ V2, suggesting MW rotation is quite flat in the vicinity R - Rq. 



3.1.1. Inner Galaxy Rotation 

Because of extinction and the relatively low velocity dispersion of the interstellar gas 
(o" v ~8 km s" 1 ), ubiquitous radio ISM gases (H 1, 12 CO) are better tracers of MW rota- 
tion in the plane than are stars, which are veritable bullets with cr v ~30 km s . Inside 
of the Solar circle, the rotation curve is directly observable by isolating the part of the 
spectral profile emitted by gas at the tangent point, where the LOS passes closest to 
the Galactic centre. At this point R = Rosin£, and the observable quantity in Eqn. Q] 
becomes the circular velocity at this radius v(Rosin£) = v t + vosin^.Success at isolating 
emission from gas there depends on 1) the method used to locate v t within a profile, and 
2) the ability to distinguish it from the innumerable clouds, expanding shells, continuum 
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and self-absorption features and other small-scale features that display their own kine- 
matics. The best metho ds model the profile near v t wi th some empirical or physically- 
motivated function (e.g. lRohrfs & Kreitschmannll 19871) rather than si mply placing v t at 



some repeatable threshold within the profile (the Tg half-maximum; iLuna et alJ l2006. 



for example). As for 2), only high-resolution data like those produced by the IGPS 
can help. The two best H i rotation curves for R < Rq on the market today are ex - 



tracted from VGPS (1 st quadrant, o r QI) and SGPS (QIV) data: iLevine et all (2008); 



McClure-Griffiths & Dickey (2007). These are c ompared in Fig. [2] after rescaling to 



our adopted constants in Sec. |2l The QIV curve of iMcClure-Griffiths & Dickevl (120071) 
illustrates the benefits of high-resolution data: it is densely sampled and makes good 
use of the high-resolution SGPS data by avoiding continuum ab sorption features. How- 
ever their curve is some ~1 km s" 1 above that of ILevine et al. (2008), a difference that 
comes from modelling the profile with a complementary error function, which re stricts 
v t to be down the wing near the half-maximum in the profile. ILevine et all d2008h use 



a model profi le which is more p hysically consistent with circular motion at the tangent 
point (that of lCelnik et aDll979h . and locates v t higher up the wing at about 90% of the 
peak, where it is generally found in synthetic observations of hydrodynamics galaxy 
simulations. They also measure QI with VGPS data, and their results show the long- 
known asymmetry between the QIV and QI curves. For their study of the z-dependence 
of v(R) they smear VGPS and SGPS data into 0?5 bins, and their curves at z =0 thus 
suffer from higher scatter than lMcClure-Griffiths & Dickevl d2007l) . The CGPS has yet 
to be used to measure QI rotation for 50° < £ <90°(R =6.4-8.4 kpc) but should provide 
an outstandingly clean curve if care is taken to avoid H i continuum- & self-absorption. 
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Figure 2. MW rotation curve between R =3-14 kpc. Inner Galaxy data are from 
the VGPS, SGPS surveys (H i terminal velocity method) and the Ou ter Galaxy f rom 
the CGPS (H ii regions). The two large grey squares are points from lXuel d2008l) . 



3.1.2. Outer Galaxy Rotation 

Outside of the Solar neighborhood, MW rotation in the disc is traceable to great dis- 
tances with the radio ISM, out to the edge of the H i disc (R ~20 kpc). Because the term 
dv/dr — > at these distances, the edge of the H i surface density as measured with Eqn. 
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[T]is very sensitive to the shape of the rotation curve (lLockmanll2002l) . In turn, the slope 
of the far outer curve measured with Eqn. Q] using discrete kinematic tracers depends 
strongly on the adopted values of constants Qo and Rq, and on uncertainties in Q and R. 
In particular, uncertainties in R for spiral arm tracers like H n regions bias part of the 
curve to appear to rise like solid-body rotation (Q ~ const lBinmy & Dehnenll 19971) : in 
other words a rising circular speed curve could be observed when it is in fact constant 
or falling through the outer Galaxy. This effect is plainly visible in Perseus and Cygnus 
arm H u regions plotted in Figure [2] appearing in each arm as a rising v(R) <x R trend. 
This bias is undo ubtedly responsible f or the rise observed in the classic Outer Galaxy 
rotation curve of iBrand & Blitzl dl993 ) traced with bright H n regions. Indeed it must 



be altering the measured shape of all rotation curves constructed principal ly from spiral 



arm t r acers. Generally , contemporary published Outer Galaxy curves (e.g. lRusseil et al. 



20071 : iHou et aLll2009h are mainly compiled from stellar d istances and veloc ities from 
across the literature, with little new data (one exception is lRusseil et al.ll2007r) . These 
curves generally show constant circular speed (i.e. Q oc R~ l ) well out to the edge of the 
Milky Way's stellar disc (~15 kpc). Near-constant c ircular spee d has been confirmed by 



modelled halo potentials from SDSS observations (IXud 120081) out to an extraordinary 
R =50 kpc where it declines to only ~200 km s" 1 . 

A feature common to all Outer Galaxy curves based on discrete kinematical trac- 
ers is a large scatter in R and (to a lesser extent) in O, which overwhelms any detail 
in the curve of the kind we see in the Inner Galaxy. To see these details, 1-arcminute 
CGPS H i and 12 CO line data can be a great advantage, measuring the systemic veloc- 
ity from the thin edges of shells that surround the ionized gas and stars. Where both 
H i and CO shell edges are seen, one gets two measures of vios which are perpen- 
dicular to any motions from the shell's expansion along the LOS. To demonstrate this 
advantage we will refer to measurements of velocities and distances of some 100 H u 
regions across much of the CGPS (90° < t <193°) (from iFoster et alJl2010L for this 
review we will refer t o these tracers often). In this sample 33 Perseus spiral arm ob- 



jects are also found in lBrand & Blitzl ( 19931) (excluding those with t =170-190° due to 



the 1/sin^ projection effect) and in that paper, these show a mean angular velocity of 
Q. ave =24.67±1.04 km s" 1 kpc -1 (1-cr). Using CGPS data two velocities each ( 12 CO 
& H i) are determined for 41 Perseus arm H n regions in the same longitude range, 
resulting in Q. ave =24. 13±0.69 km s _1 kpc -1 . In the CGPS we not only find velocities 
for 8 more objects but the scatter about a common D has been reduced by 33% Q. 

The Outer Galaxy rotation curve v(R) obtained with all 100 regions is shown in 
Fig. [2] Fitted to all data in the figure (except v t from 2007 and H n regions between 
t =170-190°, marked with open symbols) is Eqn. |2] below, which give s a reasonabl e 
approximation to observed external galaxy rotation curves (M31 for e.g. Brandt 1960): 



o3/2« 

n (R) = max nV2n (2) 

Rmax [1+2 (R/RmaxT] 

The fit is very robust and shows that linear rotation v(R) rises steadily from R =3- 
8 kpc (n =0.95), peaking at the Solar velocity v max ~ vq =255 km s" 1 at R max =9.3 kpc. 



'Additionally in lFoster et all 1 2010) distances to OB stars within the H i and l2 CO shells surrounding each 
H ii region have been recalculated using a common MK Luminosity class and a common (B - V) colour 
calibration, reducing the scatter in R by 20% 
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Thereafter it declines very gently to 235 km s" 1 atR =20 kpc. Obviously Eqn. [2]cannot 
capture details of non-circular motions in the Inner Galaxy ("wiggles") and the sharp 
discontinuity preceding Perseus Arm objects in the Outer Galaxy, almost certainly the 
signature of a shock (see Sec. I3.2I ). The slope of Eqn. |2]is undoubtedly altered by it. 

As 21 cm surveys of the MW (e.g. GALFAH i, ASKAP) evolve to higher reso- 
lutions and sensitivity so too the measured Galactic rotation curve will evolve towards 
a better sampled, lower scatter version of what we see today. As we advance into the 
next decade, VLBI parallax distances and proper motions to mor e kinematic tracer s 



will sharpen the Outer Galaxy curve and its details tremendously (Re id et alj |2009a). 



It is imperative that we now begin modelling the smooth underlying rotation curve and 
density wave perturbations separately in a multi-parameter model; not an easy task, but 
we believe the next avenue of immediate advancement in MW rotation curve studies. 
This will also help to answer the long-standing question of the asymmetry between the 
QI and QIV rotation curves. Can this be accounted for purely by spiral structure? 

3.2. Spiral Dynamics: Density Wave Motions in the MW Disc 

The choreographed movement and organization of the old stellar disc into a set of 
nested elliptical orbits creates regions of enhanced stellar density, forming a non-axisymmetric 
gravitational potential A an order of magnitude or two smaller than the total potential 
Vq. The pattern in the stars circulates as a spiral density wave (hereafter DW) through 
the disc, and the response of the ISM gas is vigorous and complex, with A{R)/Vq as 
little as 3-5% producing sharp shocks in the gas and forming many seco ndary and ter- 



tiary s piral armlets downstream of the potential minimum. As shown by Roberts et al. 



dl975D a galaxy's Hubble type and degree of development of spiral structure is deter- 
mined by the strength of shocks in its ISM, which in turn is governed by two key DW 
parameters: the pattern's angular speed of rotation Q p and the winding (pitch) angle of 
the spiral arms i. Q. p is one of the most poorly constrained parameters in all of Galactic 
astronomy. Estimates are based on either stellar kinematic tracers giving Q. p =20- 
30 km s" 1 kpc" 1 (see iGerhardl I2008L for review), or on linking periods of terrestrial 
climate change with the Sun's passage t hrough sp iral arms, giving 12-2 km s" 1 kpc" 1 
(Gies&Helsel|2005|; ISvensmarkl 2006), though lOverholt et al. 1 (120091) calls this link 



into question. It is not even clear if a single constant value for Q p can adequately de- 
scribe the Milky Way spiral. Meas urements of i are across the entire range of SBa 
to SBc galaxies (6-20° IValleel l2008h . probably in part because no distinction is made 



between what is being measured: the DW potential (stellar) arms, shocks (H i self- 
absorption & du st), and star-formation arms (H n gas and OB stars) are all expected to 
have different i ( Gittins"& Clarke! 120041) . No conclusive confrontation of observations 



with the non-linear DW (shock) theory has been presented so far, and where the MW fits 
into the spectrum of spiral galaxies remains a mystery. Are we more like galaxy M33 
(smooth continuous ISM flow; poorly defined spiral structure) or M81 (sharp shocks 
and arms in ISM and many secondary gaseous arms in the outer disc)? 

3.2.1. Density Waves in the Inner Galaxy 

We can get an initial idea of t he pattern speed in the Inner Galaxy from the recent 
GLIMPSE results (Figure 14 in IChurchwell et al.ll2009h . Star counts towards the Scu- 



tum and Centaurus arms show peaks, whereas none are seen towards the Sagittarius & 
Norma arm tangencies, suggesting that these two arms are predominantly ISM concen- 
trations. Theory predicts such secondary arms form from the gas' dynamical response 
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to a principal 2-arm stellar DW. A family of 12 hydrodynamics galaxies are built with 
combinations of Q. p =12,16,20 and i =8,10,12,14°, spanning today's range of con- 
tentious values for these parameters. Surface density images of these simulations are 
shown in online figures at www.novell.brandonu.ca/~fostert/DWModels.pdf. In them, 
the Sun is the bright dot located at Rq =7.8 kpc; its phase with respect to the pattern is 
set by the parallax distanc e to W30H {I =134°, r =2 kpc), the location of the Perseus 
arm shock (IXu et al.ll2006h . All models form a 4-arm ISM pattern beyond some radius 
which is relatively insensitive to pitch angle (R >11 kpc, 8.5 kpc & 6 kpc for Q. p =12, 
16 & 20 km s -1 kpc -1 respectively). The observation of gaseous arms deep in the Inner 
Galaxy at R ~3-5 kpc (beyond the bar) suggests then that Q p >20 km s" 1 kpc -1 there. 

The best evidence for density waves in the Inner Galaxy are the wavy irregularities 
in H i terminal velocities (~ ±10 km s -1 ) with Galactic longitude (see Fig. [3]). Here 
as well spiral structure, inextricably linked with spiral dynamics, is staring us in the 
face and crying out to be properly (quantitatively) modelled with theory in a parameter 
fitting analysis. However, Fig. [3] shows that v t is too poorly resolved, and that the 
signatures of linear and shock theories are too similar to conclusively choose between 
them with the current data. It is not clear whether more careful measurements of v t will 
unambiguously resolve the sharp cusps in v t that shock theory predicts, one of the only 
features at the tangent point distinguishing it from the linear version. Nonetheless, the 
use of high-resolution 21 cm data as illustrated by iMcClure -Grifht hs & Dickey! (120071) 



is the correct way to begin to resolve DWs in the Inner Galaxy. 

The most natural way to model DWs is by comparing theoretical 7g {(, vios ) dia- 
grams to observed ones, as theory predicts both ingredients that shape these diagrams: 
the (perturbed) density and total velocity (assuming constant spin temp erature) . Qual - 
itative compar is ons of the £ ,vj,qs diagram were made very early by iBurtonl (Il97ll) : 
Robertsl (fl972h : ISawal (Il978h and others, and once again, it was found that linear and 



shock version of DW theory produce very similar diagrams with only subtle differ- 
ences (ISawalll978h . Fig. [3] (right) shows the observed QI (,vlos diagram (Leiden H i 
data) compared to that simulated from the model galaxy with Q p =20 km s -1 kpc -1 , 
i =14°, the model we find that best fits the Leiden observations in a simple x 1 - 

z(T° B bs - T'™ d f minimization. 

3.2.2. Density Waves in the Outer Galaxy 

In the Outer Galaxy (QII & QUI) the velocity anomaly associated with Perseus spiral 
arm H n regions has been long known, and is seen in published rotation curves today 
between R =9-11 kpc as a cluster of objects underneath the smooth v(R) ( s ee Fig . 
|2] The explanation in terms of the shock was proposed early on by iRobertsI dl972h . 
but the idea was never really explored further from there, probably because the scatter 
in rotation curves based on H u regions prevents a positive identification of a sharp 
feature like a shock. However, the Outer Galaxy rotation curve traced with CGPS 
H ii regions clearly shows this velocity discontinuity between the local & Perseus arm 
objects at R =8.9 kpc (see Fig. |4j left). To compare this to, Q.(R) is simulated from 
our family of synthetic galaxies by picking any point a known distance from the Sun's 
location and assuming its entire LOS velocity projection is described by Eqn. Q] The 
results of picking thousands of these points show the type of rotation curve one might 
measure for such a galaxy. Fig. |4] (left panel) is the simulated rotation curve from 
model Q p =16 km s -1 kpc -1 , i =12° indeed showing that a shock would appear as a 
discontinuity in Q. extending under the basic rotation curve at R =8.9 kpc, just as is 
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90 70 60 50 40 30 20 l^l^^H^flli^l M^^^^^H 

Longitude ^(degrees) 80° 60° 40° 80° 60° 40' 



Figure 3. (Left) Inner Galaxy (QI) terminal velocities v, predicted from circular 
rotation (dashed), linea r DW theory (dotted ) and non-linear theory (solid), compared 
to measured ones from lLevine et al.l ([2008) (circles). (Right) Synthetic T\, ((, vws) 
diagram (right) compared to the Leiden H i survey (left) for i =22-90°. Both panels 
are calculated from the same galaxy model Q. p =20 km s _1 kpc 1 and i =14°. 



seen in the Perseus arm H n regions. We have found that only simulations with Q. p =16 
place the shock exactly at their location; Q. p = 12 models place the shock well ahead of 
the H ii regions, and Q. p =20 models do not form shocks at all beyond the Solar radius. 

We distinguish a group of more distant H ii regions at a constant Q =20 km s" 1 kpc" 
which trace a second optical arm some ~2 kpc beyond Pe rseus. The presence of sec- 
ond op tical arm beyond the Perse us arm was proposed bylKimeswenger &~W einberger 
(Il989u . and stars in it traced by iNegueruela & Marcol d2003l) who refer to it as the 
"Cygnus arm". In our sample some 20 objects are distinguished from the Perseus arm 
members by not only their distinctly larger distances (R ave ~H kpc), but by 1) their 
relatively high latitudes (mean 2?5) compared to the Perseus arm, 2) their correspon- 
dence with an H i concentration similarly higher up in the distant warped midplane (see 
Fig. [6]>, and 3) their appearance in a larger "kinematic ring" (Q. - Q.o) oc |vLos/sin 
In our simulated rotation curve Cygnus arm objects appear behind a secondary shock 
beyond the Perseus arm (Fig. [4]). The presence of both star-forming arms at R > Ro 
has implications for the pattern speed Q p , if we are to assume that the H ii regions 
that trace them form from shocks. The right panel shows a birds-eye view of model 
Qp =12 km s" 1 kpc" 1 , i =12°. This model generates multiple strong secondary shocks 
behind the primary one in the Perseus arm (at R =8.9 kpc) that explains the presence 
of H ii regions in the Cygnus and Outer arms (yellow and green respectively). The 
Q. p =16 models do not generate shocks beyond the Perseus arm shock, and Q. p =20 
models show no strong shocks form at all beyond R = Ro, making the presence of any 
massive star-formation regions outside the Solar circle not as easily explained. 

It is likely then that the M W has a radially falling spiral pattern speed of the kind 
seen in some external spirals ( Meid tlitai] 120091) . with Q. p >20 from the bar's edge to 
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Figure 4. (Left) Outer Galaxy Q,(R) observed with 100 H n regions in the 
CGPS, plotted with the simulated rotation curve observed in a synthetic galaxy 
(model Q p =16 km s _1 kpc 1 , i =12°). (Right) Surface density map of model 
Q p =12 k m s _1 kp c" 1 , ; =12°with 100 CGPS H n regions overlaid not fitted. See 
http://www.novell.brandonu.ca/~fostert/DWModels.pdf for larger figures. 



R ~ Rq, dropping to Q. p ~16 at R ~9 kpc and to Q. p ~12 km s -1 kpc -1 at R >11 kpc. 
Admittedly this conclusion is drawn from somewhat simple comparisons made for the 
sake of this review, but we believe demonstrate what is possible to discover in con- 
fronting DW theory with kinematics of the Outer Galaxy. Large-scale comparisons 
of this kind with new observations have today fallen behind the ever-increasing num- 
ber of theoretical simulations. The power in new high-resolution 21 cm surveys will 
be their ability to resolve the smooth intercloud medium between the clouds, allowing 
H i cubes that show a less-confused view of large-scale Galactic mo tions like rotation, 
density wave streaming and turbulence to be created. iRohlfs! ( 19741) recognized this 
and applied a minimum- brightness boxcar typ e filter to isolate the large-scale H i "stra- 
tum" from the "clouds". iLevine et all (l2006al) use an unsharp mask on LAB H i data 
to enhance the visibility of spiral structure in surface density maps. A similar more ad- 
vanced filter based on the background separation technique in Kot hes & Kertonl (|2002|) 
is demonstrated with CGPS data in Fig. [5] With any filter it is important to use high 
angular resolution observations in order to resolve the narrow gaps between clouds. 
Creating new expressions of the neutral ISM images will ultimately resolve the subtle 
differences between density wave theories. 



3.3. Rolling Motions 

Spiral arms in most directions in the MW exhibit strange velocity gr adients perpendic- 
ular to the plane. These were first observed in the 2 1 cm line by Oort (1 19621) and coined 
"rolling motions". The term comes from the arms' tilted appearance in H i latitude- 
velocity slices Tb (b,VLos)\ at positive latitudes the 7# centroid of the arm appears at 
higher negative velocities, whereas down through the equator and below (b < 0) the 
centroid smoothly shift s to more positiv e velocities. Rolling is a lso observed in the 
molec ular ISM ( 12 CO,[Wouterlooj| 1981 ). and in local young stars ( Brosche & Schwar] 
1984). Discrete objects like H n regi ons also "roll", exh ibiting the same overall gradient 
across the plane as the H i arms do ( Foster et alj|2010l ). and are therefore likely tracing 
the same motions. Typical gradients are dv/db ~2-4 km s" 1 per degree of latitude, or 




Figure 5. A 10° wide Tb (i, vlos ) plot from CGPS H i data separated into "cloud" 
(left) and "intercloud medium" (right) maps using a spatial filter. Very little coherent 
large-scale structure remains in the cloud component, in contrast to the intercloud 
map at right which is readily compared to DW models. 



about dv/dz ~20 km s 1 per kpc in the Perseus arm. Rolling motions have not yet been 
obse rved in external galaxie s. 

I Yuan & Wallace! (1 19731) first described rolling as apparent "pseudo-motions" cre- 
ated by our top-down view of the arms caused by their tilt towards up (due to the warp 
of the midplane), and differential circular rotation. For example, in QII the upward 
warp causes a given Tb contour of an arm at higher latitudes to be observed at a greater 
distance (and thus at a more negative velocity) than the same contour at lower latitudes. 
This geometric effect is su btle, and is reproducable in 3 D models of the Galactic disk 
that include the warp (e.g. iFoster & MacWill iams!l200r3) . however, a modern flat rota- 
tion curve and figure for the warp produces sma ller apparent gradients (<1 k m s _1 per 
degree) than observed. This was also shown by iFeitzinger & Spickerl d 19851) . opening 
the possibility that rolling motions are partly physically real. They (1986) propose tum- 
bling pockets of H i gas falling from high in Galactic Fountains, themselves caused by 
active star formation in t he arm's middle. 

The lYuan& Wallace] model seems to be supported by observations in the Outer 
Galaxy. If rolling is entirely due to observed differences in circular motion (includ- 
ing rotation) with latitude, then gradients in the arms should disappear in the cardinal 
direction i =180°. This is observed for both the Perseus & Cygnus H i arms (see 
Fig. ©. However, the latitudes b of Perseus arm H ii regions are constant (b ave =0?1) 
with increasing distance R, yet these objects also "roll" with the H i arm. Rolling in 
Perseus must therefore be entirely physical and not a viewing effect. The magnitude of 
the rolling gradient in Cygnus arm objects is the same as in Perseus, so rolling in the 
Cygnus arm could also have essentially nothing to do with its displacement above the 
midplane. However H n regions associated with this arm do show a relationship of R 
with latitude b, just as one would expect if the warp were involved! 

As the reader can see there are more inconclusions than conclusions. How much 
of the observed gradients in the MW are real, and how much apparent? A good 3D MW 
model that includes the warp fitted to the observations will be necessary to separate real 
from apparent rolling motions, after which the fountain theory can be further tested 
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VELOCITY (km/s) [LSR] Longitude I (degrees) 



Figure 6. (Left) Tg (b, vlos) slice from the DRAO 26-metre H i survey 

(Higgs & Tapping 2000), showing the Perseus (v L0S 50 km s _1 , b -0°), the 

Cygnus (-70 km s~\ +3°) and the Outer H i arm (-90 km s" 1 , +2.8°). The triangle 
marks the b, vios of Sh2-141, seen within the H i concentration that iFoster et al.1 
(2010) identify with the Cygnus arm. (Right) Gradients measured across QII for 
both arms. Functions dv(i)/db = -2.44sin(f)+0.03 and -2.35sin(f)+0.06 fitted to 
the Perseus (solid line) and Cygnus arms (dotted line) respectively are drawn. 



or other explanations sought. Are rolling motions unique to the MW, and therefore 
likely illusions caused by our viewing point? Future 21 cm observations of external 
galaxies will tell us. The EVLA is the only instrument with resolution <1 km s" 1 and 
~l"on the sky sufficient to scrutinize external spiral arms with 50 pc spatial resolution 
at D=10 Mpc. Nearby galaxies in the THINGS survey (see this volume) would be good 
starters, and in the future ASKAP will be able to extend our observations. 



4. Galactic Structure 

An executive summary of the MW might go something like this: the MW is a barred 
spiral galaxy with a short boxy bulge in its centre, and a stellar bar about twice as long. 
The old stars form a principally 2-arm DW pattern in the Inner Galaxy, giving rise to at 
least 4 ISM arms that are active in star formation from R ~4-14 kpc. The 4-arm pattern 
winds through most of the Inner and all through the Outer Galaxy with the same pitch 
angles of i ~12°, and the arms are connected by minor branches and spurs throughout, 
including one that the Sun currently resides in. Inside the Solar circle the H i disc is 
remarkably flat, while outside the disc is twisted nearly along the Sun-centre line into 
a gentle warp above the midplane in one hemisphere and below it in the other. 

4.1. Milky Way Spiral Structure: Methods and Results 

Mapping the spiral structure of the Milky Way historically is still an exciting challenge 
that continues to this day. The historic 21 cm map by Kerr-Westerhout is based on 1) 
velocimetric deconvolution of the observed Tb {£, vr.ns) brightness dist ribution using 



circular rotation and a constant spin temperature. iLevine et all (|20 06a) recreate only 



the Outer Galaxy this way using the LAB H i data, and a minor refinement to the veloc- 
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ity field. Although a bit better than the historic map, the same distorted picture emerges 
with pitch angles i =21-25°, differing wildly from those more directly measured with 
kinematic tracers. These maps demonstrate that H i is a poor tracer of Galactic spi- 
ral structure when used like this. Deconvolutions based on the iDame et al.1 ( 200lh CO 
survey are cleare r (e.g. iNakanishi & Sofuell2006l . find a more reasonable i =12-15°). 
Poh l et all d2008) use a non-circular velocity field based on a hydrodynamics (HD) 
galaxy simulation, itself fitted to the observed NIR luminosity distribution. The result- 
ing Inner Galaxy map is more detailed than any previous one, showing a distinct bar in 
the centre angled to the Sun-centre line by <p =20°, and two arms emanating from the 
bar's ends at R ~3.5 kpc and winding outward by i ~1 1?5. Thus far no study has applied 
such a model to deconv olution of the o bserved H i distribution, though the approach is 
promising as shown by iGomezI (120061) in reconstruction of observed simulations. The 
more flexible approach of 2) modelling the obser ved Tr (I, b, v) distribution is demon- 
strated in Fig. |7|(right). An empirical DW model dFoster & MacWil liamsll2006l) created 
the synthetic £ — v plot at bottom, and 12 adjustable parameters were varied in a proper 
least-squares fit to the data (top). Only two DW arms are modelled. Tight winding an- 
gles of i =6°are found, primarily because the model treats the Cygnus and Outer arms 
as one. The computational cost of this multi-parame ter model-fit ting approach is not 
prohibitive as it is for more realistic MHD models (Gomez 2006). Empirical models 
can also include components that have been difficult to incorporate in smoothed particle 
hydrodynamics models, such as the warp, magnetic fields, an d rolling motions. 

The paradigmatic 4-arm map of iGeorgelin & Georgelinl d 19761) is based on 3) di- 
rect spatial mapping using 268 H ii regions with photo metric distances and/or kine- 
matic distances based on RRL and Ha velocities. The IGeorgelin & Georgelinl map 
has achieved some fame (cited at least 550 times ) and some infamy. It is the starting 
point for the also-f amous electron density map by lTaylor & Cordesl ( 1993 ). and the dust 
e mission model by iDrimmel & Spergel (1200 ll) . Compared to modern maps the pattern 
in IGeorgelin & Georgelinl is unusually (and worringly) clean, and has prompted many 
to propose the MW is among the 10% of spirals with a so-called "grand design". While 
it is generally accepted that there is indeed a coherent spiral pattern to be mapped, 
H ii regions and giant molecular clo uds in the MW sport a truly messy pattern (c.f. 
Paladini et al.ll2004l : iHou et al.ll2009 ). revealed by "photoshopping" out the fitted arms 



drawn on published H u region maps. One wonders if H ii regions will ever reveal a 
"g rand design" spiral st r ucture in the MW , if one exists at all. Recent H ii region maps 
by iRusseil et al.l d2007l) : IHou et all ( 20091) fit 2-4 logarithmic or variable pitch angle 
spiral segments, with 4 arms at i =10-13°always fitting better than two with tighter 
i =6-8° windings. Of course, any model with twice the free parameters should fit better, 
but 4-arm models also reproduce known arm tangency directions. Within the likely un- 
certainties of the fits (±20-30%; see below), 4-arm models show the same pitch angle 
arm to arm and negligible variation of i with radius (this however is not clear as few 
papers offer proper uncertainties on any fitted spiral arm parameters!). A least-squares 
fit to CGPS H ii regions weighted by uncertainties in both coordinates InR & cf> yields 
i =12.0±2.6°and 12.6±4.1°for the two major star-forming arms in QII (see Fig. [7). 

Another method, 4) fitting models to Galactic features observed on the sky (e.g. 
arm tangency directions) yields information on the number of arms as well as large 
symmetric structures like the disc scale length , thick ness and flaring, and the warp. An 
example is the model of IDrimmel & Spergell (l200ll) based on C OBE observations of 
NIR (stars) and FIR (warm dust) emission. IDrimmel & Spergell intriguingly find that 
an m =2 arm pattern dominates the nonaxisymmetric NIR emission. This is supported 




Figure 7. (Left) Spiral structure in QII and part of QUI from CGPS H n re- 
gions. Fitted pitch angles are: Perseus arm (red diamonds) i =12.0+2.6°; Cygnus 
arm (yellow triangles) 12.6±4.1°, and Outer arm (gr een squares) 5.6+3.9°. ( Right) 
Tb {(, vlos ) plot at b =0°from the 26 m H i survey b y lHiggs & Tapping] (|2000) com- 
pared to the model of iFoster & MacWilliamsl ([2006) fitted to it quantitatively. 



by Inner Galaxy star counts from the GLIMPSE survey, which show no enhancement in 
star n umbers towards the ex pected tangency of the Sagittarius arm ~45°< t <51°(Fig. 
14 in lChurchwell et alj|2009l) . A visualization of the MW constraine d with GLIMPSE 
H i & CO Galactic plane surveys, and other observations (Fig. 15 IChurchwell et al. 



2009) shows the Perseus and Scutum-Centaurus arms are the principal stellar arms, 



while Sagittarius and Norma are two star-poor ISM arms. The Norma arm wraps around 
from its origin at R =4 kpc to bec ome the "Outer arm" in QII & QUI, and then is 
identified with the distant H i arm of iMcClure-Griffiths et all (120041) in QUI & QIV. 

C omparison of the GLIMPSE map with the recent synthesis of MW spiral struc- 
ture by Vallee (2008), and with CGPS H n regions, shows some fundamental differ- 
ences. Vallier (l2008l) draw the Sagittarius and Norma arms attached to the bar's ends, 
whereas the Scutum-Crux and Perseus arms first appear outside the zone of the bar. In 
the Outer Galaxy the w inding angle and distance to the Perseus arm agree extremely 
well in all three maps. Vallee (2008) show the Cy gnus arm as the Outer Galaxy ex- 
tension of the Norma arm, however, it appears muc h farther from the Sun (6-7 kpc 
between I =130-150°) than H ii regions and stars (e.g. Neguer uela & Marcoll2003L an d 
Fig. |7J this work) show it to be (4-5 kpc). It is possible that its distance in IValled is 
kinematically estimated and overly large due to the presence of a second shock beyond 
the Perseus arm in QII. The Cygnus arm is missing altogether from the GLIMPSE map, 
but would appear halfway between the Perseus and Outer arms on there. Given the per- 
vasive confusion in the literature between it and the Perseus arm, this is not surprising. 
MW spiral structure is then far from a closed subject, and as pointed out by 



Churchwell et al.1 is deserving of much more attention from the community. In particu- 



lar the relationship between the Inner and Outer Galaxy patterns needs to be clarified. 
How do the stellar and ISM arms from GLIMPSE relate to the Outer Galaxy pattern? 
As computational speeds increase, new attempts at determining MW spiral structure 
in H i by modelling and comparisons with DW theory, and deconvolution with self- 
consistent dynamical models may be profitable. As mor e VLBI parallax distances can 
be measured towards fainter masers dReid et al.l l2009ah . the precise pattern (or lack 
thereof!) traced out by star-forming regions will become clearer in the decade to come. 
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4.2. The Warp 

The pioneering 21 cm observations of the MW disc in the 1950's showed the remark- 
able flatness of the Inner Galaxy disc, and in dramatic contrast, the warp of the Outer 
H i layer. A warped midplane is a commonality shared by most galaxies whose H i 
extends past their o ptical discs. Our neighbors M31 and M33 both exhibit warps (e.g. 
Che min et al. I2009h . The warp remains the least understood large-scale component of 



MW structure. It begins at R max ~ 1 .5/?o inside the stellar disc. The Sun-centre line also 
happens to nearly bisect the warped halves of the H i disc, with the half in the North- 
ern sky much more pronounced in amplitude. In QI & QII, H i profiles in directions 
€ < 120° and b >0 are characteristically shaped by the LOS passing t hrough, above, and 
then p lunging back through the flared outer disc at oblique angles dFoster & Routledge 
2003b . A good analytic figure for the warp is written as sinusoidal in azimuth (f> and 



linear in R (m =1), with a small quadratic (m =2) term: 



z(R, <p) = ai(R - R max )sm(f) + a 2 (R - 7W) 2 sin 2 (3) 



Binnev & Merrifield ( 1998) give a\ =1/6, a 2 =1/120 & R max =11 kpc, appropriate for 
a flat rotation curve with Rq =8.5 kpc & Q.q =27.5 km s" 1 kpc -1 . The observed z > 
displacement of the midplane is maximum towards i =80° , whe r e the H i layer extends 
more than 20 degrees (3 kpc in z) from the equator. iLevine et al.1 d2006b) fit m =0, 1 & 2 



harmonic modes with the same equation, finding the m =1 mode begins inside the Solar 
circle and dominates to R =16 kpc, where the m =2 mode begins to appear, introducing 
the asymmetry in the North and South halves of the warp. Observations that 3 out 
of 4 warped galaxies have no visible companions suggest that warps are not tidally 
excited by satellites. However, when dark matter is included in simulations with small 
companions (<2% total mass), the "wake" left by the orbiting com panions in the DM 
halo in duces vibration modes in the disc. In the MW, simulations by lWeinberg & Blita 
(2006) show a warp with 3 modes (m =0 1 & 2) can be excited by the Magellanic clouds 



passing through the DM halo. 

The CGPS shows other essentially unstudied modes in the disc's topography with 
wavelengths of 5-20 kpc spectacularly well. These are collectively referred to as "scal- 
loping", and give the Outer disc a crinkly, corrugated appearance. 
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